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ABSTRACT 

The quiet Sun photospheric plasma has a variety of magnetic field strengths going from zero to 1800 G. 
The empirical characterization of these field strengths requires a probability density function (PDF), i.e., a 
function P(B) describing the fraction of quiet Sun occupied by each field strength B. We show how to combine 
magnetic field strength measurements based on the Zeeman effect and the Hanle effect to estimate an unbiased 
P(B). The application of the method to real observations renders a set of possible PDFs, which outline the 
general characteristics of the quiet Sun magnetic fields. Their most probable field strength differs from zero. 
The magnetic energy density is a significant fraction of the kinetic energy of the granular motions at the base 
of the photosphere (larger than 15% or larger than 2 x 10 3 erg cm" 3 ). The unsigned flux density (or mean 
magnetic field strength) has to be between 130 G and 190 G. A significant part of the unsigned flux (between 
10 % and 50 %) and of the magnetic energy (between 45 % and 85 %) are provided by the field strengths larger 
than 500 G which, however, occupy only a small fraction of the surface (between 1% and 10%). The fraction 
of kG fields in the quiet Sun is even smaller, but they are important for a number of reasons. The kG fields 
still trace a significant fraction of the total magnetic energy, they reach the high photosphere, and they appear 
in unpolarized light images. The quiet Sun photosphere has far more unsigned magnetic flux and magnetic 
energy than the active regions and the network all together. 

Subject headings: Sun: fundamental parameters - Sun: magnetic fields - Sun: photosphere 



1. INTRODUCTION 

The interior of the supergranulation cells does not seem 
to be mag netic according to the stan dard full disk magne- 
tograms. iLivingston & Harvevl (119751) and ISmithsonl (119751) 
improved the standard sensitivity to discovered weak mag- 
netic signals in the cell interiors. These elusive signals are 
known as Inter-Network fields (IN), Intra-Network fields or, 
simply, quiet Sun fields. Characterizing these quiet Sun 
fields is important since they may play a role i n the physics 
respons ible f or the g lobal solar magnetism j UnnoJ [19591 
Stenfl oiri98l lYi et all 1 19931 Islnch ez Alm eidalll998l 120041 
Schriiver & Title 2003). The quiet Sun could easily contain 
most of the photospheric unsigned magnetic flux and mag- 
netic energy, simply because it occupies most of the solar sur- 
face ( more than 90 %, even at solar maximum: lHarvev- Angle! 

Our understanding of the quiet Sun magnetic fields 
has improved during the last years thanks to advances 
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tions (e.g. ICattaned . 

IStein & Nordlundl 120021 IVoglerl 120031 IVogler et ail 120051) . 
The new spectropolatimeters make it possible to detect 
the low polarimetric signals characteristic of the IN fields. 
Sophisticated diagnostic techniques are needed to interpret 
those signals, which depend in a non-trivial way on the 
(a priori unknown) physical properties of the magnetized 
plasma. For example, opposite polarities often coexist in 
the resolution elements giving rise to very asymmetric line 
profil es requiring special methods of diagnostic (see, e.g., 
Sanchez Almeida & Lite! 120001 iLitesI 120021) . The numer- 
ical simulations provide guidance for the magnetic field 
topology to be expected, a knowledge critical for a proper 
interpretation of the observables. 

The magnetic field strength is one of the key parame- 
ters to characterize the magnetism of the quiet Sun. Its 
study has a rathe r long tradit i on, with a period of unsuc- 
cessful attempts (lUnnol Il959t iHoward & B hatnagar Il969t 
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Stenflol (11982) and Fauro bert-Scholll ( 119931) based on Hanle 
depolarization measurements. The observed Hanle signals 
agree with a volume-fillin g mixed-polari t y magnetic field 
of con stant field strength ijS tenfld Il982l iFau robert-S cholil 
1 19931 IFaurob ert-Scholl et al.1 ll995tnFaurobert et al.l 1200 ll 
ISterm^^t^niT997ir The value of the estimated mean field 
strength has increased upon refineme nt of the diagnostic 
techni ques. The most recent e stimat e bv lTruiillo Bueno et alJ 
(2004) and Bommi er et al.l (120051) are consistent with a 
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constant field of 60 G (or a mean field of 130 G assuming 
an exponential distribution of magnetic field strengths). 
Most measurements of the quiet Sun magnetic fields are 
based on interpreting line polarization signals caused by the 
more familiar Zeeman effect. They show a rather different 
magnetism, with structures occupying a small fraction of 
the solar surface and having intense fields in the range of 
the hectogauss (hG) and the kilogauss (kG). Actually, the 
field strength depends in a systematic way on the Zeeman 
splitting of the spectral line employed to measure. The works 
using visible lines deduce t he pre s ence of kG magnetic fields 
iGrossmann-Doerth et alJ 119961: ISanchez A lmeida <^LjJgs| 
2000; Socas-Navarro & Sanchez Almeida 2002; 

Domf nguez Cerdena et alJ l2003albl) . while measurements 
based on infrared (IR ) lines infer sub-kG fie ld strengths 
JLin & Rimmeld Il999t iKhomenko et alJ 120031) . This ap- 
parent inconsistency among the various measurements 
of field strength, plus the guidance of numerical simula- 
tions of magneto-convection, suggest that the quiet Sun 
has a continuous distribution of field strength s going all 
the w ay fr om G to 2 kG (ISanchez A lmeida & Lites| 
2000: ISocas-Navarro & Sanchez AlmeidaT 120031 
Truiillo Bueno et al. 2004). Therefore, each one of the 
individual measurements mentioned above is biased, pro- 
viding only a partial view of the magnetic fields existing in 
the quiet Sun. We think that this biased information can be 
cleaned up and assembled to provide a reasonable estimate of 
the intrinsic magnetic field strength distribution. This work 
describes a method to recover such information. 

Due to the variety of field strengths, they are best charac- 
terized by a probability density function PDF. The magnetic 
field strength PDF, P(B), is defined as the probability that a 
point of the atmosphere chosen at random has a magnetic field 
strength B. Alternatively, it describes the fraction of quiet Sun 
occupied by fields with a strength B. The PDF approach to 
characterizing the quiet Sun field strength is both convenient 
and powerful. First, one can directly compare observations 
with numerical simulations, which predict P(B). Second, it 
concentrates basic information on the quiet Sun magnetism, 
e.g., the first moment (B) is connected with the unsigned mag- 
netic flux density 1 , whereas the second moment provides the 
magnetic energy density (B 2 )/(87r), 

/>oo 

(B)= / BP(B)dB, (1) 
Jo 

(B 2 }= / B 2 P(B)dB. (2) 
Jo 

There are several observational magnetic field strengths PDFs 
in the literature (see the references cited in the previous para- 
graph). As we pointed out above, all of them are known to be 
severely biased. On the one hand, those based on Hanle sig- 
nals miss the strong fields, since the Hanle depolarization does 
not vary for field strengths above a threshold that is typically 
small e r than a few hundred G (see, e.g.,|L andi DeglTnnocenti 
119921 IStenflol U99H iTruiillo Buenol 120011) . On the other 
hand, the Zeeman measurements show only a small fraction 

1 Strictly speaking, the unsigned magnetic flux density across a plane is 
given by the average of the unsigned component of the magnetic field in the 
direction perpendicular to the plane. Consequently, we are not precise when 
using the term unsigned magnetic flux density for (B) , which is the average 
of the magnetic field vector modulus. There is no possible misunderstanding, 
though. The term is employed consistently throughout the text, and the actual 
unsigned flux is never mentioned except in this footnote. 



of the existing magnetic stru ctures (see, e.g., IStenflolfl982t 
ISanchez Almeida et alJl2003h . The Zeeman signals tend to 
cancel when opposite polarities coexist in the resolution ele- 
ments. We detect only the residuals left when the cancellation 
is not perfect (e.g., when the two polarities do not have exactly 
the same magnetic flux, when strong Doppler shifts separate 
the absorption produced by the two polarities, when the tem- 
perature imbalance between polarities creates net signals, and 
so on). Here we try to get rid of all these observational bi- 
ases to produce an empirical P(B). We approach the problem 
by assuming a flexible yet reasonable shape for the expected 
PDF. This guess is based on recent numerical simulations of 
magneto-convection, as well as on the latest Hanle and Zee- 
man observations. The actual P(B) is set by forcing the PDF 
to meet several observational and theoretical constraints. So 
far as we are aware of, this is the first quiet Sun PDF satis- 
fying magnetic field strength observations going all the way 
from G to 2 kG. (The p remises were already put forward by 
ISanchez Al meida 2004, though.) The success of the trial is 
only partial since rather than a single PDF we derive a fam- 
ily compatible with the present observations. However, the 
approach is shown to be viable, so that future more precise 
observations will allow us to narrow down the range of possi- 
bilities. 

The paper is organized as follows; § |2] describes how to 
piece together all the ingredients needed to estimate an unbi- 
ased PDF. Some of the more tedious details are completed in 
Appendix [A] The properties of the PDFs satisfying the ob- 
servational constraints are extracted and analyzed in § [3] In 
particular, we compare the semi-empirical PDFs with theo- 
retical PDFs coming from numerical simulations of magneto- 
convection. Finally, § 0] summarizes the main conclusions, 
discusses their consequences, and points out the natural way 
to proceed. 

2. SYNTHESIS OF THE PDF 

In order to facilitate the reading of this detailed section, we 
start off by listing the thread of the argument. 

1 . The unbiased quiet Sun PDF can be approximated by 
a linear combination of the biased PDFs i nferr ed from 
Zeeman signals and from Hanle signals (§ 12. U . 

2. The shape of the Zeeman PDF is obtained from obser- 
vations. The shape of the Hanle PDF is not constrained 
by observations, and so we take it from numerical sim- 
ulations of magneto-convection (§ 12. 1> . 

3. The Hanle and Zeeman signals come from different 
photospheric heights. One needs a transformation to 
express them in a common height, in other words, a 
recipe for the varia tion with height in the atmosphere 
of the PDFs (§l2~2l. 

4. Piecing together this information, we work out a semi- 
analytical PDF depending on three free parameters. 
O ne ne eds only three independent observables to set it 

(§E3- 

2.1. The shape of the PDF 

The Zeeman induced polarization changes the sign depend- 
ing on the orientation of the vector magnetic field with re- 
spect to the observer. When magnetic fields with varied di- 
rections coexist in the resolution elements, the net polariza- 
tion tends to cancel. The lack of polarization is interpreted 
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FIG. 1. — PDFs inferred from Zeeman signals of Fe I A6301.5 A, 
Fe I A6302.5 A, Fe 1 A15648 A, and Fe I A15652 A observed simulta- 
neously in a quiet Sun region. The two PDFs come from the same data 
set. The solid line represents a MISMA inversion i Domrneuez Cerdena et al. 
2005 ) whereas the dashed line corresponds to a Milne-Eddington inversion 
i Sanch ez Almeida et allUffiS) . The MISMA PDF is shown some 100 km 
above the base of the photosphere, at a height corresponding to the Milne- 
Eddington inversion. The PDFs are normalized to the fraction of quiet Sun 
producing Zeeman polarization signals (~ 1.5 %). 



(non-local thermodynamic equilibrium). Hanle based stud- 
ies are therefore restricted to lines formed high in the atmo- 
sphere, where the collisions are low enough for the scatter- 
ing events to be frequent. However, we need Hanle signals 
formed as low as possible, so that they can be compared with 
Zeeman signals produced in the low-mid photosphere. The 
Hanl e signals of Sr I A4607 A a re formed particularly low 
(e.g. iFaurobert-Scholl et alJl!995l) . In addition, this line al- 
lows a simple treatment of the usuall y very complex syn- 
thesis of the Hanle signals (see, Faurobert et al.l l200lL and 
also the discussion in § I2.3> . Sr I A4607 A is chosen to con- 
strain the PD F for these two unusual p roperties. As it is ar- 
gued in, e.g.. ISanchez Almeidal (120 05). the Hanle signals of 
Sr I A4607 A do not constrain the shape of the PDF. This 
limitation is overcome here by parameterizing the shape of 
Ph(B) using an analytic function. The Hanle signals have been 
interpreted with approximations using Dirac (5-functions, 
exponen tial fu nctions, and Max wellia n distributions (e.g., 
Stenfl ol 119821 iFaurobert-SchollI Il993t iTruiillo Bueno et alJ 
2004; Sanchez Almeida 2004), but these shapes do not fit in 
our needs. First, we expect a continuous distribution of field 
strengths, which discards Dirac (5-functions. Second, the ex- 
ponential fall off with field strength, 



as absence of magnetic fields, so that the Zeeman measure- 
ments underestimate the volume occupied by magnetic fields. 
This bias increases with decreasing field strength, since the 
Zeeman signals scale with the field strength. The Hanle in- 
duced depolarization does not suffer the cancellation of the 
Zeeman signals. However, the Hanle induced depolarization 
is insensitive to magneti c field strengths larger than a few 
hundred G (e g.. lLandi DeglTnnocenta|1992t lStenflolll994t 
Truiill o Buenoll2001l ISanchez Almeidall2005l) . Thus the bi- 
ases of the Hanle and Zeeman measurements turn out to be 
complementary. Keeping in mind this complementarity, one 
can approximate the true unbiased PDF, P{B), as a linear com- 
bination of the biased PDF inferred from Zeeman signals, 
Pz(B), plus the biased PDF that fits the Hanle signals, P H (B), 



P{B) = W P H (B) + (\- W )P z (B), 



(3) 



with Ph(B) and Pz(B) normalized to unity, and the weight 
w ~ 1. A full derivation of this equation is given in Ap- 
pendixlAl 

Regarding Pz(B), we take the shape of the biased PDF 
measured by Dommguez Cerdena et al. (2005). It is based 
on a realistic semi-empirical modeling of 2280 Stokes pro- 
files of Fe I A6301.5 A, Fe I A6302.5 A, Fe I A15648 A, 
and Fe I A 15652 A observed simultaneously in a quiet 
Sun region. These four spectral lines are complemen- 
tary since those in the visible spectral range sense pref- 
erentially the kG fields whereas the IR lines tend to 
detect hG field strengths (ISanchez Almeid a & Lites| l2000t 
ISocas-Navarro & Sanchez Alm eida 20021 120031) . The empir- 
ical Zeeman PDF shape adopted here is not very different 
from the one derived by mea ns of Mil ne-Eddington fits to 
the same Stokes profiles (see ISanchez Almeida et alJ 120031 
Dommguez Cerdena 2004J). The two PDFs are shown in Fig- 
ure \l\ Both are rather broad with a local maximum at kG 
fields. The differences are unimpo rtant for the present work, 
and they are discussed in detail by Dommguez Cerden a et al.l 
(2005). 

The Hanle effect modifies the polarization produced in scat- 
tering processes, and so it only affects lines formed in NLTE 



PH(B)=—rexp(-B/(B)) : 



(4) 



is too pronounced according to the PDFs resulting from nu- 
merical simulations of magneto-convection (see, e.g., the dot- 
ted line in Figure^). The Maxwellian distribution can be dis- 
carded with the same argument, since the fall off is even more 
pronounced than the exponential one. Inspired by the numer- 
ical simulations of magneto- convection develop ed with the 
MuRAM code l lVoglerll2003l IVogler et"alll20 05). we choose 
a log-normal distribution Ph(B), 



Ph(B) = 



1 



expf-anfi-lnBo) 2 /^ 2 ]: (5) 



with Bo and a two parameters related to the two first moments 
of the distribution, 



and 



(B)=B exp(a 2 /4), 



(B 2 )=B 2 sxp(a 2 ). 



(6) 



(7) 



Figure [2] compares numerical PDFs from IVogler] £2003) 
with fits using log-normal functions and exponential func- 
tions. (All the fits were performed for B < 500 G.) The two 
plots correspond to an initial vertical field of 10 G (Figure^) 
and 50 G (Figure and the parameters of the log-normal 
PDFs are a = 2.5, B = 5 G, and a = 1.7, B = 38 G, respec- 
tively. The log-normal fits are much better than the exponen- 
tials, and we choose them to represent Ph(B). A comment is 
in order. The log-normals tend to zero when B — >0 G. Actu- 
ally, they reach a maximum at B = B max with, 



B max = B () exp(-a 2 /2). 



(8) 



Counterintuitive as it sounds, the property of having a max- 
imum is to be expected for the true quiet Sun PDF. When 
the field strength is small enough, the magnetic fields are 
easily modified and randomized by the convective motions. 
Those weak fields are expected to have an isotropic distribu- 
tion of field directions, with the three components of the mag- 
netic field being independent. Suppose that the magnetic field 
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FIG. 2. — Fits of quiet Sun PDFs from n umerical simulations of magneto- 
convection (the solid lines; Voaler 2003) using log-normal functions (the 
dashed lines) and exponential functions (the dotted lines). The log-normals 
are far superior. The upper and lower panels differ in the initial field strength 
of the simulation; 10 G (a) and 50 G (b), respectively. The PDFs correspond 
to the base of the photosphere. 



strength of these weak fields has a typical value B, different 
from zero but otherwise arbitrary. For a point of the atmo- 
sphere to have B = the three components of the magnetic 
field vector must be zero simultaneously. This event is highly 
improbable since the three components are independent. Con- 
sequently, it is to be expected that P(B) — >0 when B « B,. 
The drop at small fields is also clearly seen in some numerical 
simulations of magneto-convect ion (e.g., Fig ureEb; Figure 3 
in lCattaneoll999t Figure 5.25 in Vogl erl2003l) . 

2.2. Vertical variation of the PDF 

The photosphere is strongly stratified, with the gas pres- 
sure decreasing exponentially with the height in the atmo- 
sphere. Part of the pressure that maintains the mechanical 
balance of a magnetic structure is provided by the magnetic 
pressure, which has to decrease to keep up with the pressure 
drop. Since the magnetic flux is conserved, the decrease of 
field strength comes together with an expansion of the field 
lines and an increase of the volume occupied by the magnetic 
fields (see, e.g., Spruit 1981a). These two effects modify the 
PDF. As we pointed out above, the Hanle signals are formed 
high in the photosphere. In order to use them to constrain a 
PDF at the base of the photosphere, one needs to know how 
the PDF varies with height. Working out such variation is a 
non trivial problem, closely connected with the so-called ex- 
trapolation of the photospheric fields to the chromosphere and 



(a) view from a side 




FIG. 3. — (a) Cartoon with two magnetic fiuxtubes. The field lines of one of 
them reach the mid photosphere (z = 250 km). The field lines of the second 
one bend over and return below this photospheric height, (b) A view from 
the top of the first fluxtube. The thick solid circles outline the sections of the 
fluxtube at z = km (small circle) and z, = 250 km (large circle). Note that 
for the fluxtube to expand a number of other concentrations around it cannot 
reach z = 250 km. The section of those that must return are represented by 
dotted circles in the cartoon. 



corona (e.g.. lSakuraJT98^lSchriiver & Titldl2 003>. Here we 
develop a crude approximation, which suffices for the kind of 
exploratory investigation carried out in the paper but which, 
eventually, will require further refinement. 

We start off by introducing a new symbol for the PDF at a 
height in the atmosphere z, p(B,z), so that the PDF discussed 
so far corresponds to the base of the photosphere 2 , 

P(B)=p(B,0). (9) 

Our estimate is base d in the so-called thin fluxtube approxi- 
mation (e.g.. lSpruiil981albl) . In this approximation, the mag- 
netic field strength of an isolated static structure surrounded 
by field free plasma drops with the height in the atmosphere 
as the square root of the field free gas pressure. In our case, 

B(z)~B(0)exp(-z/H B ), (10) 

with magnetic field strength scale height Hb twice the gas 
pressure scale height, and so, independent of the field 
strength. The cross-section of one of these thin fiuxtubes S 
increases with height as the field strength decreases. The need 
to conserve the magnetic flux imposes, 

S(B(z),z) =S(B(0),0) exp(z/H B ), (11) 

so that the product B S does not depend on z- The expansion is 
made at the expense of the field free plasma surrounding the 
magnetic concentration. In our case there is no such field free 
plasma, so that an increase of the volume occupied by one 
magnetic structure requires shrinking the volume occupied by 
others (see the two magnetic concentrations in Figure ^Jt). 
Which physical parameter determines why some concentra- 
tions expand whereas others shrink? The magnetic forces 
responsible for the expansion or contraction scale with the 
square of the magnetic field strength. Then it is to be expected 
that the expansion of the strongest fields prevails, making it 

2 The height z = corresponds to the base of the photosphe re in the one - 
dimensional mean quiet Sun model atmospheres (e.g. Maltb vet alJH98fl) . 
where the continuum optical depth r c = land the total pressure is ~ 1.2 x 
10 5 dyn cm . 
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difficult for the weak fields to reach high photospheric lay- 
ers. We use this idea as an ansatz to proceed. For a magnetic 
concentration to expand from S(B(0),0) to S(B(z),z), its field 
strength has to be the largest in the area S(B(z),z)) around it 
(see Figure^). If in average all magnetic concentrations have 
the same section, then a concentration will survive if and only 
if it has the largest field strength among the n — 1 concentra- 
tions around it, with 



S(B(z),z) 



(12) 



S(B(0),0)' 

see Figure [3J). The probability that one concentration has a 
field strength smaller than B at z = is, 



M(B)= / P(B')dB'. 



(13) 



Since the field strengths of different concentrations are inde- 
pendent, the probability that n - 1 have fields smaller than 
B is M(B)"~ l and, consequently, a magnetic field strength B 
reaches z with a probability, 



F P (B) = M{Bf 



(14) 



This probability of survival Fp(B) allows us to write down the 
PDF at z, namely, 



p(B(z),z)dB(,z) = 



S(B(z),z) 
S(B(0),0) 



F P (B(0))p(B(0),0)dB(0), 



(15) 

where we have taken into account that the atmosphere oc- 
cupied for those magnetic fields that survive increases as 
s(b(z),z)/S(B(0),0). Using equations ©, dT0|, O, and 
111 21 . one can rewrite equation d!5l > in a compact way, 

p(B,z)=n 2 P(nB)[ P{B')dB'}"~ X . (16) 
Jo 

It is not difficult to show that fjo(B,z) is properly normalized 
even for non-integer n, which we interpret as a sign of self 
consistency in the derivation of equation d!6i . 

Equation J16I provides the PDF at any height given n and 
the PDF at height equals zero. We apply such recipe to the 
Hanle signals of Sr I A4607 A when they are formed as close 
as possible to the disk center (i.e., as low as possible in the 
photosphere). T his lowest height of fo rmation turns out to 
be about 250 km (Faurobert-Scholl etafll 19951) . On the other 
hand, a H B — 230 km is obtained from the stratification of 
B in th e inversions performed by iDommguez Cerdena et alJ 
(2005), which renders, 

n = exp(250km/// B )~3. (17) 

Figure0]shows the variation with height predicted by equa- 
tion dl6> for an exponential P(B) (a), and a log-normal P(B) 
(b). One can compare the original PDFs (the dashed lines) 
with the PDFs at z = 250 km (the solid lines). Note how the 
maximum of P(B) is shifted to stronger fields as one moves 
up in the photosphere. This is to be expected since the proba- 
bility of survival of weak fields is very low. 

2.3. Physical constraints 

Onc e Pz(B) has been set according to the prescription in 
§ 12.11 The analytical PDF in equation Q depends on three 
free parameters; the weight w, and the two parameters a and 
Bq defining the shape of the Hanle PDF (equation 0). They 
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FIG. 4. — Examples of the change with height in the photosphere of an 
exponential PDF (a) and a log-normal PDF (b). The plots in the first column 
show the probability of survival Fp(B). The plots in the second column repre- 
sent the PDF at z = (the dashed lines), and the PDF at z=250 km (the solid 
lines). (The PDFs correspond to those shown in Figure Eb.) 



can be set using three independent observables. Unfortu- 
nately, only one observable plus one lower limit are available. 
First, the fraction of the photosphere producing the observed 
Zeeman signals is a few per cent ( e.g.lLin & Rimmeielll99 
Sanchez Almeida & Lites 2000; Dommguez Cerdena et al 
2003a; Lites & Socas-Navarro 2004). This fraction represents 
a lower limit to the true percentage of qu iet Su n with strong 
fields, which constrains 1 - w (equation IA17I ). In particu- 
lar, the signals used to estimat e Pz(B) arise from some 1.5% 
of the surface (Dommguez Cerden aet alJl2005T) . Taking this 
value as our lower limit, 1 — w > 1.5 x 10~ 2 , and so, 



w< 0.985. 



(18) 



Second, the Hanle depolarization signals of Sr I A4607 A, 
Q/Qo, depend on the PDF approximately as, 



Q/Qo = / $*B,z)W B (B)dB, 
Jo 

with the kernel Wb{B) given by, 



W B (B) = l-~ 



1 2 H 



47l 



5Vl +7 | l+4 7 2 
B 

7h = 



(19) 



(20) 



50 G' 



see. lFaurobert et afl ( l2001l) . ISanchez Alm eida (2005), and ref- 
erences therein. Using equations d!6l > and d!7i . equation d!9l 
can be rewritten as, 



3B 



Q/Q ~9 W B (B)P(3B)[ P{B')dB'} dB. (21) 
Jo Jo 

The work by Truiill oBueno et aT] ( 120041) gives our best esti- 
mate of Q/Qo, 

Q/Q = 0.41 ±0.04, (22) 
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with the error bars representing the standard deviation among 
different observations made at a fixed heliocentric angle. Ful- 
filling equations i2li and 1221 warrants producing the ob- 
served Sr I A4607 A Hanle depolarization. We will refer to 
them as the Hanle constraint. 

Lacking of additional observables to determine a single 
P(B), we resort to a general consideration of theoretical na- 
ture. The magnetic energy associated with the quiet Sun mag- 
netic fields has to be smaller than the kinetic energy of the 
granular motions, otherwise the non-magnetic numerical sim- 
ulations of solar granulation would not b e so successful in 
reproducing all kinds of observables (see iStein & Nordlundl 
H998tlAsplund et al.l2000l) . In other words, 

X«l, (23) 

if x is the ratio between the magnetic energy and the kinetic 
energy of the granular motions at the base of the photosphere, 

{ ^ 1 = x\{pu 2 ). (24) 

Using the typical convective velocities u and densities p at 
z = k m (p ~ 3 x 10" 7 g cm" 3 , u ~ 3 km s -1 ; Stein & Nordlund 
119981 Figure 5), one gets a kinetic energy density of some 
1.3 x 10 4 erg cm" 3 . This value allows us to rewrite equa- 
tion J24> as, 

X =(fl 2 )/(580G) 2 . (25) 

3. RESULTS 

The PDF in equation Q d epends on three free parameters. 
As it is discussed in § 12.31 we do not have three indepen- 
dent constraints to determine a single P(B). There is one gen- 
uine constraint (equations 1211 and 1221 ) and two lower limits 
(equations 1181 and 1231 ). This problem rules out selecting 
one PDF for the quiet Sun magnetic field strength. Instead, 
we are forced to explore the set of PDFs compatible with the 
observables, and then to distill properties common to all of 
them. The strategy of the search and the common properties 
are described next. 

Given a, we produce a grid of 18 x 30 P(B) with 0.83 < 
w < 1, and 5 G < B < 150 G. It suffices to include all the 
possibilities allowed by the observational and theoretical con- 
strains. For each PDF we compute the total magnetic energy 
in units of the kinetic energy (\ in equation 1251 ). and the 
Sr I A4607 A Hanle signals to be expected (equation 1211 ). 
The magnetic energy and the Hanle signals as a function of 
w and Bo are shown in the plots of the first column of Fig- 
ure|5](column a). The total energy is represented by the solid 
contours. The dashed lines trace those pairs w - Bq compat- 
ible with the observed depolarization i22\ . The shaded strips 
around these dashed lines correspond to the error bars of the 
depolarization. The five plots portray four different values of 
a (0.5, 1, 1.75, and 2.5), plus the case of exponential Hanle 
PDFs (equation |4]). The first column of Figure [5] shows that 
the magnetic energy is always important (say, larger than 15 % 
of the kinetic energy of the granular motions). This conclu- 
sion turns out to be independent of Bo and a provided that w 
fulfills the limit <QJQ . 

The third column in Figure[5]shows specific PDFs for three 
pairs w - Bq producing the observed Hanle signals. The 
dashed line, the dotted line and the solid line correspond to 
the circle, the cross and the square in the first column, re- 
spectively. They differ mainly because of w, a parameter con- 
trolling the importance of the tail of strong fields, with little 



influence on th e co re of the PDF at low field strengths. As 
we discuss in § 12.11 the PDFs present a maximum in this core 
of small B. (The most probable magnetic field strength dif- 
fers from zero.) Since the Hanle constraint determines Bo 
almost independently of w (the dashed lines in Figures |5^ 
are nearly horizontal), the position of this maximum is con- 
trolled by a (see equation (8)). Figure [6] shows how the 
magnetic field strength of the maximum B max increases with 
decreasing a to become B max ~ 180 G when a — > 0. This 
limit can be explained keeping in mind that P(B) tends to a 
Dirac <5-function 5(B — Bq) when a — > (see equation Q). 
A Dirac <5-f unction at B ~ 60 G accounts for the Hanle con- 
straint (l22l (ITru iillo Bueno et aT1 l2004t iBommier et alJl2005l 
Sanchez Almeida 2005). According to the recipe that we use 
for the variation with height of the PDFs (equations 1161 and 
El), a S(B-60 G) at the height where Sr I A4607 A is formed 
corresponds to 5(B- 180 G) at z = 0, which explains the limit 
fl max ~ 180 G. 

In order to illustrate some other details of the PDFs fulfill- 
ing the observational constraints, we select a representative 
one having a few additional desirable properties. It will be 
denoted as the reference PDF. In principle, there are no ar- 
guments to prefer a specific B max , and a wide range of val- 
ues is compatible with the observables (Figure |6j. How- 
ever, we tend to think that small Z? max PDFs are more real- 
istic, since they seem to be prefe rred by the numerical sim- 
ulations of magne to-convection ( Emonet & Cattaneo: 1200 It 
IStein & Nordlundl l2002t IVo glerl 120031: IVogler et alJ 120051) 
Small B max implies large er, but a cannot be arbitrarily large 
since the magnetic energy of the PDF increases with a (Fig- 
ure |5^), and the magnetic energy must comply with the 
limit d23l >. We choose a = 1.75, which yields B max ~ 13 G. 
This selection together with w = 0.94 leads to the reference 
PDF shown in Figure It has an energy 28% of the kinetic 
energy (\ = 0.28), and an unsigned flux density (B) ~150 G. 
The dashed line in Figure^ also shows the reference PDF at 
250 km, computed according to the recipe in equation J! 6b . 
The latter has an unsigned flux density of 120 G, similar to 
the flu x of the exponential PDF used bv ITruiillo Bueno et alJ 
(2004) to account for the Hanle constraint in equation J22l> . 
The magnetic energy at z = 250 km is three times smaller than 
the energy at the base of the photosphere. FigurelTJ) compares 
our reference PDF and simulated PDFs from Vogler (2003) 
with initial fields of 10 G (the dotted line), 50 G (the dashed 
line), and 200 G (the dotted-dashed line). These 3D simu- 
lations solve the MHD equations for a compressible partly 
ionized plasma. After a purely hydro-dynamic transient, an 
initial homogeneous and vertical field is introduced. The 
processing and amplification of these seed fields lead to ad- 
ditional magnetic flux, which amounts to 25 G, 100 G and 
270 G, for the three different initial fields. The reference PDF 
turns out to be a mixture between the 200 G PDF (B < 600 G) 
and the 50 G PDF (1000 G < B < 1500 G). Note that the 
two numerical PDFs have a small broad local maximum at 
large magnetic fields (B ~ 1500 G). Our semi-empirical PDF 
also has a maximum at large fields, except that it is more pro- 
nounced, and it occurs at higher fields (B ~ 1700 G). The 
existence of such maximum at hight field strengths can be 
understood as the result of a magnetic intensification process 
from weak to strong fields. The magnetic energy piles up at 
the largest possible field strength, set by the value whose mag- 
netic pressure exceeds the gas pressure of a non-magnetic at- 
mosphere in mechanical equilibrium (~ 1800 G at z = km). 
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FIG. 5. — (a) Magnetic energy density of the PDFs as a function of the free parameters w and Bo ( me solid contours). The labels give the magnetic energy in 
percent of the kinetic energy of the granular motions. The dashed lines show those pairs w - Bq fulfilling the Hanle constraint, with the shaded strips representing 
to the range of observed depolarizations. The different plots portray different values of the parameters a, plus the case of an exponential Hanle PDF included for 
reference (bottom; in this case Bo = (B)). The three symbols correspond to the three symbols in the second column, and to the three PDFs of the third column 
(see below), (b) Fraction of quiet Sun unsigned magnetic flux density (the dotted line) and of magnetic energy density (the solid line) in the form of strong fields 
(B > 500 G). They are represented versus their filling factor a(B > 500 G). The curves correspond to pairs w - B fulfilling the Hanle constraint (the dashed 
lines of the first column), (c) PDFs corresponding to the three symbols in the first and the second columns. The solid line corresponds to the square (w = 0.9), the 
dotted line to the times symbol (w = 0.98), and the dashed line to the bullet (w = 1.0). 
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FIG. 6. — Magnetic field strength at which P(B) is maximum versus a. 
The different curves have different w; 1, 0.98 and 0.9, for the solid line, the 
(barely visible) dotted line, and the dashed line, respectively. 
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FIG. 7.— Reference PDF with a = 1.75 and w = 0.94. (a) PDF at z = km 
(the solid line) and at z. = 250 km (the dashed line), (b) Comparison of P(B) 
at z = (the solid line) with PDFs from numerical simulations of magneto- 
convection by Voalei 1 2003) having initial field strengths of 10 G (the dotted 
line), 50 G (the dashed line) and 200 G (the dotted-dashed line). 



(A detailed explanation of the mechanism can be found in 
iDomfnguez Cerdena et 31120051 1 

We are particularly interested in understanding which part 
of the quiet Sun magnetism is provided by strong fields which, 



among other things, are responsible for most of the observed 
Zeeman signals (see AppendixlEl. and can be studied in unpo- 
larized light (see below). In order to quantify the contribution, 
we use the first moments of the PDF. Specifically, the fraction 
of quiet Sun occupied by strong fields, 

a(B>B*)= / P{B)dB 1 (26) 

the fraction of unsigned magnetic flux, 

0(B>B*)=^- f BP(B)dB, (27) 
{B) J B * 

and fraction of magnetic energy, 

1 f°° 

e(B >B*) = —j B 2 P(B)dB. (28) 

Choosing a somewhat arbitrary threshold between weak and 
strong fields 3 , we use 

B* = 500 G, (29) 

for the plots in the second column (b) of Figure|5] They show 
the fraction of unsigned flux (the dotted line) and magnetic en- 
ergy (the solid line) versus the volume occupied by the strong 
fields. Only PDFs producing the observed Hanle depolariza- 
tion are considered (those corresponding to the dashed lines 
in the first column of Figure|5J. The dots, crosses and squares 
also correspond to the dots, crosses and squares in the first 
column of Figure|5] In all cases with w < 0.98 (indicated by 
the cross symbols), the fraction of magnetic energy in strong 
fields is larger than 50% (between 45% and 85% when 0.5 < 
a < 2.5), although they fill only a small fraction of the sur- 
face (between 1% and 10%). The contribution is dispropor- 
tionated. As far as the unsigned flux is concerned, the strong 
field still provide a sizeable fraction (between 10% and 50% 
for 0.5 < a < 2.5). The reference PDF (a = 1.75, w = 0.94) 
has 42 % of the u nsigned flux in B > 500 G, a figure s imilar 
to the estimates bylS^o^g^javgrro & Sanchez Almeida! (12003! 
38%) and lSanchez Almeida! (120041 36%1. 

According to a physical mechanism originally proposed 
by Sp ruitl dl976l) . the very strong fields are expected to be 
bright in unpolarized images. Part of the pressure required 
to maintain the plasma in mechanical balance is provided 
by the magnetic pressure. When the magnetic pressure of a 
particular structure is comparable to the mean gas pressure 
of the atmosphere, the structure must have low mass den- 
sity and so it becomes transparent. We see deeper through 
the magnetized structure, and deeper often means hotter and 
brighter. The mean pressure at z = is equivalent to a mag- 
netic pressure provided by a field strength of some 1800 G. 
In order to estimate which fraction of the quiet Sun mag- 
netic fields could be bright, we consider the strongest mag- 
netic fields in the PDFs, say, those larger than two-thirds of 
the maximum value, or larger than 1200 G. Figure [8^ shows 
a(B > B*) versus B* for three representative PDFs fulfill- 
ing the Hanle constraint. The reference PDF corresponds 
to the solid line. The filling factor in bright features is al- 
ways very small, e.g., a(B > 1200 G) =2.5% in the refer- 
ence PDF. Yet, these fields are responsible for 25 % of the 
unsigned flux and 50 % of the magnetic energy of the quiet 

3 This figure is not far from the value of a magnetic field strength yielding 
equipartition betw een magnetic energy density and kinetic energy density; 
see equation 1251 with \ = 1 . 
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FIG. 8. — Filling factor (a), fraction of unsigned magnetic flux (b), and 
fraction of magnetic energy (c) in field strengths stronger than B* . The three 
curves correspond to three PDFs fulfilling the Hanle constraint with a = 1.75 
and w = 0.94 (the solid line, and also the reference PDF), a = 1.75 and w = 1 
(the dotted line), and a = 1 and w = 0.94 (the dashed line). 



Sun magnetic fields. ISanchez Almeida et al 1 d200l detected 
(G-band) bright points in the quiet Sun covering some 0.7 % 
of the supergranular cell interiors. Although the percentage is 
smaller than the surface coverage in strong fields mentioned 
above, the two figur es may not necessarily di sagree. As it is 
pointed out by Sanchez Almeida et al. (2004), the observed 
value is strongly biased and it represents a lower limit. Only 
the largest magnetic co ncentrations pop up ab ove the dark in- 
tergranule background (Title & Berser 1996) and, moreover, 
not all kG are expe cted to be bright (Sanchez Almei d~et al.l 
l207villVoglei<20na . 

So far the value of Q/Qo given in equation (I22> has been 
used as the Hanle constraint. In order to explore the de- 
pendence of the results on Q/Qo, we also tried Q/Qq ~ 
0.6, which is the Hanl e depolarization inferred and used by 
iFauro bert etafl ( E00l . The main effect is decreasing the 
value of Bo needed to reproduce the observations, and so, de- 
creasing the magnetic field where P(B) is maximum (B max , 
see equation (8)). This shift of B max can be achieved without 
a large variation of the unsigned magnetic flux and magnetic 
energy. This is due to the fact that the tail of hG and kG 
field strengths is the mayor contributor to the flux and energy, 
and this tail is not very sensitive to Bo, and so, to B max . For 
instance, the reference PDF (w = 0.94, a = 1.75) has B max = 
3 G, (B) =70 G, and \ - 0-17 when Q/Q ~ 0.6, whereas 
Bmax = 13 G, (B) =150 G, and X - 0.28 when Q/Q ~ 0.41. 
This moderate change in (B) and \ must be compared with 
the case where the tail of strong fields is not present (w = 1 in 
equation Q). Then (B) oc Z? max and \ ^ ^m ax (equations ]6), 
(7). and JHJ), which implies a very large decrease of the un- 



signed flux (a factor of four) and energy (a factor of twenty) 
to fit a depolarization signal of Q/Qo — 0.6. 

The unsigned magnetic flux and magnetic energy assigned 
by our semi-empirical PDFs to the quiet Sun is far larger than 
the unsigned flux and energy in the form of active regions and 
the network, even at solar maximum. The reference PDF has 
(B) ~ 150 G and (fi 2 ) 1 / 2 ~ 305 G. Consider the regular full 
disk magnetograms, which are mostly sensitive to the active 
regions and the network (see §Q. At solar maximum, these 
magnetograms show a total unsigned magnet ic flux across the 
solar surface of some 7 x 10 23 Mx (e.g.. lHarvev-Anglell993L 
Chapter 12). Dividing this magnetic flux by the area of the 
solar surface, the unsigned flux density in the form of active 
regions and the network (B)ar turns out to be, 

~ 12 G ~ 0.08 (B), (30) 

where we have used for (B) the unsigned flux of the reference 
PDF. If all these active region magnetic fields are intrinsically 
strong with B — 1500 G, then 

(B 2 )ar * 1500 G (B) AR ~ (135 G) 2 ~ 0.2 (B 2 ). (31) 

4. CONCLUSIONS 

The quiet Sun magnetic fields may be central to understand 
the global magnetic properties of the Sun, but they are not 
well characterized from an empirical point of view. They 
present a variety of physical properties and so their descrip- 
tion must be carried out in terms of probability density func- 
tions (PDFs). This work is focused on characterizing the PDF 
of the field strength, which describes the fraction of quiet Sun 
occupied by each field strength. We use the symbol P(B) for 
the PDF at the base of the photosphere 4 . The PDFs existing 
in the literature provide only a partial view - each individual 
measurement samples only a small range of field strengths 
(see §Q. Our paper proposes a semi-empirical P{B) that tries 
to get rid of all the biases and so, for the first time, it cov- 
ers all the range of possible field strengths from B = G to 
B ~ 1800 G. It combines Hanle effect based measurements, 
Zeeman effect based measurements, and some general ansatzs 
inspired by numerical simulation of magneto-convection and 
the thin fluxtube approximation (§ |2). The procedure is of 
general nature and eventually, it will allow us to determine a 
single PDF characteristic of the quiet Sun magnetic fields. So 
far the observational and theoretical constraints only provide 
a range of PDFs. The constraints used in this work are: (a) 
reproducing the Hanle depolarization signals of Sr I A4607 A, 
(b) having a filling factor for strong fields larger than the fill- 
ing factor inferred from the observed Zeeman signals, and (c) 
having a magnetic energy density smaller than the kinetic en- 
ergy density of the granular motions. 

The properties common to all P(B) fulfilling these con- 
straints are analyzed in §|5] They provide a qualitative picture 
of the quiet Sun magnetic field plus some quantitative results. 

- The magnetic energy density is a significant fraction 
of the kinetic energy density of the granular motions. 
Larger than some 15% of the kinetic energy at the base 
of the photosphere. 

- Most of the surface is occupied by weak magnetic 
fields. We find that magnetic fields weaker than 500 G 
fill between 90% and 99% of the surface. 

4 Where the continuum optical depth at A5000 A is one, and the total 
pressure ~ 1.2 X 10 5 dyn cm -2 . 
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- The most probable field strength is not zero. This result 
has a novel consequence. The existence and position 
of such maximum is a direct indication of the presence 
of large amounts of unsigned magnetic flux in the quiet 
Sun. Its presence represents a specific prediction of the 
work, and it may be tested observing Hanle signals of 
lines with assorted sensitivities to weak magnetic fields. 

- The magnetic energy tends to be concentrated in strong 
fields independently of the parameters chosen for the 
PDF. Between 45 % and 85 % of the energy is in the 
form of fields stronger then 500 G which occupy only 
a small fraction of the quiet Sun (between 1 % and 
10%) . The unsigned magnetic flux of these strong 
fields represent between 10% and 50% of the quiet 
Sun unsigned flux. Since the strong fields reach the 
u pper photosphere more easily thant the weak ones 
(§ 12.2b . they may be particularly influential in connec- 
tion with the c hromospheric and coronal magnetism 
(e.g. lSchriiver & Titlel2003l 

- The tail of very strong kG fields of the PDFs is expected 
to be bright in hig h spatial resolution unpq l arize d im- 
ages (Spruj3[l976[). ISanchez Almeida et ail (120041) de- 
tected (G-band) bright points in the quiet Sun super- 
granulation cell interiors covering some 0.7 % of the 
surface. The filling factors predicted here are larger 
but not necessarily incompatible. The observed bright 
points are only a fraction of the existing ones (the 
largest and brightest ones). Despite the minute frac- 
tion of the surface that it covers, this tail of kG fields 
contains a significant part of the quiet Sun magnetic en- 
ergy. They have the practical advantage with respect to 
the weak fields of showing up in unpolarized light. 

- The unsigned magnetic flux and magnetic energy as- 
signed by our semi-empirical PDFs to the quiet Sun are 
far larger than the unsigned flux and energy in the form 
of active regions and the network, even at solar maxi- 
mum. 

The sensitivity of the current Zeeman measurements allows 
us to detect spatially resolved unifor m magnetic fields of one 
G or even smal ler (e.g.. iKeller et alJll994t iLites et al]ll99rt 
iLin & Rimmele1ll999t iKhomenko et alJl2003l) . Most of the 
solar surface is occupied by magnetic fields with strengths 
larger than this sensitivity limit which, however, do not show 
up in the measurements (c.f. the Zeeman inferred PDFs in 
Figure \l\ with the reference PDF in Figure Q. The Zeeman 
signals seem to be dominated by the strong magnetic fields, 
despite the fact they only fill a small part of the solar surface. 
As we explain in Appendix 151 the reason is twofold. First, 
the Zeeman signals are sensitive to the magnetic flux, and the 
unsigned magnetic flux in the form of weak and strong fields 
is comparable. Second, all Zeeman signals tend to cancel, but 
the cancellation of opposite polarity low field strength signals 
is more effective. Typically, there are many more weak mag- 



netic field features per resolution element than strong mag- 
netic field features. Then the weak fields cancel more often 
than the strong fields, and the latter are detected in a larger 
proportion. (See AppendixlElfor quantitative arguments.) 

Our conclusions are in the vein of previous works assign- 
ing large amounts of unsig ned magnetic flux and magnetic 
energy to the quiet Sun ( ISanchez Almeida & Litesl 120001: 
Dommguez Cerdena etalj l2003aF Hjnche^Almek ii^Uri] 
20031 ISanchez Almeidal I2003t ISanchez Almeidal 12004 
Truiill o Bueno et alJ 12004 iBommier et alJ 120051) ! On the 
other hand, they seem to contradict another line of research 
indicating that the quie t Sun is not so ma g netically active 
JFaurobert-SchoIil fl993t iWang et al.l [T9951 iFaurobert et ail 
1200 It IKhomenko et alJ 120051) . As we discuss in § □ all 
previous estimates are known to be biased, and so, they 
require model-dependent corrections to go from the actual 
measurements to the solar magnetic fluxes and energies. Pos- 
sibly a significant part of the discrepancy can be pinned down 
to such corrections, and it will be cured upon refinement of 
the techniques and improvement of the statistical significance 
of the results. 

A quiet Sun as magnetic as the one advocated here may 
pose a problem of consistency with existin g hydro-dynami c 
numerical simulations of convection (e.g., Spruit et al. 1990). 
The simulations reproduc e with realism a large number of ob - 
servables (e.g.. lStein & Nordlundll998tlAsplund et al.1 2000). 
even though they do not include magnetic fields. Could the 
quiet Sun have large amounts of magnetic flux and energy 
without modifying the convective pattern in a significant way? 
In other words, what are the minimum unsigned flux and en- 
ergy required to alter the convective pattern? These are open 
questions to be addressed with realistic numerical simulations 
of magneto-convection. 

On final comment is order. We simplify the problem of 
characterizing the quiet Sun magnetic field strengths using 
a single representative PDF. However, it is to be expected 
that different constituents of the quiet Sun have very differ- 
ent PDFs. For example, we do not distinguish between the 
PDFs of gra nules and intergranule s although both numerical 
simulations dCatta neoll999t|Vogler| 2003|) and measure ments 
( Socas-Navarro etal] 12004 iTr uiillo Bu eno et all 12004) indi- 
cate a clear preference of the strong fields for the intergran- 
ules. The separation between the PDFs in granules and in- 
tergranules remains pendi ng o f future work. Similarly, the 
Zeeman PDF adopted in 32. H is taken from a specific obser- 
vation, and it may depend on the spatial resolution and other 
specificities of the particular observation. Studying such de- 
pendence remains pending too. 
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APPENDIX 

THE UNBIASED PDF IS A LINEAR COMBINATION OF BIASED PDFS COMING FROM ZEEMAN AND 

HANLE MEASUREMENTS 

The spectra used to derive magnetic field strengths from the Zeeman polarization signals have a finite spatial resolution. 
Therefore, two polarities often co-exist in the resolution elements, which cancels part of the polarization signals. As a result, 
the Zeeman measurements provide only a fraction f(B) of the quiet Sun magnetic fields with strength B. The complementary 
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1 —f(B) is interpreted as field free. Using the symbols P(B) for the true PDF, Pz(B) for the biased PDF inferred from Zeeman 
signals, and 5(B) for a Dirac (5-function, 

Pz(B) = P(B)f(B) + [l-(f)]S(B), (Al) 
with (/) the fraction of magnetized quiet Sun according to the biased Zeeman measurements, 



(/>= / 
Jo 



f(B)P(B)dB, (A2) 



and Pz(B) normalized to unity. We ignore the details of how the fraction f(B) depends on the field strength B. However, it is to 
be expected that the bias of the Zeeman PD F increases toward low field s trengths because th e polarization signals weaken with 
decreasing field strength (e.g. lUnnoll956tlLandi DeglTnnocenti & Landi Degrinnocentill973l) . In addition, the bias f(B) should 
not vary with the field strength for magnetic fields stronger than a certain limit B\ . We conjecture that given a spatial resolution, 
the cancellation of polarization signals becomes independent B when B > B\. One can summarize these general conditions to be 
satisfied by f(B) as, 

fW^if/ 1 1>b\. ^ 

Likewise, the PDF derived from Hanle signals, Ph(B), also differs from the true PDF. The difference is parametrized by g(B), 

P H (B)=^P(B). (A4) 

(«?) 

The function g(B) is also unknown, except that it has to allow Ph(B) to be a good approximation of P(B) at small field strengths. 
In addition, the Hanle signals are insensitive to field strengths larger than the so-called Hanle saturation B2 (see §[3- This fact 
allows us to assume g(B) ~ for B > B2 since no information on the true PDF is provided by Ph(B), 

1 B^O, 



««HJ B>B 2 . ^ 



For the sake of convenience, we consider Pz(B) without the (5-function. After normalization, the process renders a new Zeeman 
PDF, 

P Z (B) = [Pz(B)- (1 - (/)) 5(B)] (/)"' = j^-P(B). (A6) 



In principle, one could infer the value of P(B) from both the Zeeman PDF, 

>_(/>_ 

/OB) 



P(B) = P z (B)^l-, (A7) 



and the Hanle PDF, 



8(B) 

In fact, any linear combination of equations (IA7> and (IA8> also renders the true PDF, 



P(B) = P H (B)^-. (A8) 



P(B) = h(B)Pz(B)^L + [l-h(B)]P H (B)^- (A9) 
f(B) g(B) 



with a free weighting function h(B). If one chooses h(B) to be, 

= ' /i 



h(B)=^P-, (A10) 



then, 

P(B)=^P z (B)+[l-^-]^Lp H (B), (All) 
/1 f\ 8(B) 

or adding and subtracting (g)P H (B), 

P(B)=Vlp z (B) + (g)P H (B)+ [l-i^l-g(B)]P(B). (A12) 
/l /l 

The third term of this general expression is negligible for both small and large B, i.e., 

\-f(B)/f 1 -g(B)«\. (A13) 

When B is large then g(B ) ~ and f(B ) ~ f\ (see equations IA3I and IA5I ). The same happens for small B since f(B) / f\ << 1 
and g(B) ~ 1 (equations IA3I and IA5I ). Continuity arguments suggest that the term also vanishes for intermediate strengths, so 
that 

(f) 

P(B) ~ ^Pz(B) + (g)P H (B). (A14) 
/1 
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Using the normalization condition, 



DC 



P(B)dB = P z (B)dB = P H (B)dB=\, (A15) 

JQ JO 

the expression (IA14I can be rewritten as, 

P(B)~wP H (B) + (l-w)P z (B), (A16) 

with w = 1 - (/) //i . Note that, 

w<l-</), (A17) 

since f\ < 1 . Then the weight w has t o be smaller than one minus the filling factor that the Zeeman signals assign to the magnetic 
fields of the quiet Sun (equation IA2I ). This constraint is used in the main text. 

In short, with very unspecific assumptions on the biases affecting the observed Zeeman and Hanle PDFs, one can approximate 
the true PDF as a linear combination of the observed biased Zeeman and Hanle PDFs. 

WHY UBIQUITOUS WEAK FIELDS DO NOT PRODUCE SIGNIFICANT ZEEMAN SIGNALS. 

The sensitivity of the current Zeeman measurements would allo w us to detect a magnetic field of one G or even s maller 
if it is uniform and fi lls completely the resolution element (e.g., Kell er et aDll994t iLites et a"flll996t iLin & Rimmeldll999t 
iKhomenko et al.ll2.Q03h . According to the arguments in the main text, most of the solar surface is occupied by magnetic fields 
with strengths larger than this sensitivity limit which, however, do not show up in the Zeeman measurements (c.f. the measured 
Zeeman PDFs in Figure [0 with the unbiased PDF in Figure 0. Actually, the Zeeman signals seem to be dominated by the 
polarization created by strong magnetic fields, despite the fact they only fill a small fraction of the solar surface. Signals from 
opposite polarities existing in the resolution element cancel out, an effect invoked to explain the difference between the real and 
the observed PDF (Appendix |A). Why does this bias affect preferentially the weak fields? Why do the small magnetic field 
strength signals cancel more efficiently than those from large field strengths? 

The observed bias can be understood if the size of the unresolved magnetic structures contributing to the Zeeman signals is 
uncorrected with the magnetic field strength. (The typical size of a unipolar structure does not depend very much on its field 
strength.) Then, typically, each resolution element contains far more structures having weak fields than structures with strong 
fields. The cancellation is more effective when there are many structures that cancel signals, which explains the bias. This 
qualitative argument can be quantified as follows. A unipolar patch produces a circular polarization Stokes V signal which, to 
first order, scales with the field strength 5 . Then the Stokes V signal to be expected in a finite resolution element due to magnetic 
fields of strength B is 

V = V B^B)[N + (B)-N-(B)], (Bl) 

with N + (B) and N~(B) the number of structures with positive and negative polarities in the resolution element. The symbol Vb 
stands for an irrelevant constant, whereas /i(B) represents the typical value of the cosine of the magnetic field inclination with 
respect to the vertical direction. (We consider observations made at the disk center.) The number N + (B) can be regarded as a 
random variable, in the sense that it has different values in different resolution elements. Since it is a number of counts, the 
random variable N + (B) follows a Poisson distribution (e.g., Mart ini 19711) . whose expected value and variance are identical, 

E{N + (B)} = V{N + (B)}. (B2) 

Moreover, 

P(B) = E{N + (B)+N-{B)}(l/L) 2 , (B3) 

with L and / the size of the resolution element and size of one of the magnetic elements, respectively. (The factor [L/l] 2 is the 
number of structures required to cover a resolution element.) There is no preferred polarity, consequently, 

E{N-(B)} = E{N + (B)}. (B4) 

Equations (IB3> and JB4> lead to 

E{N + (B)} = ^-(L/l) 2 . (B5) 



If the random variables N + (B) and N (B) are independent, one can readily co mput e the mean sign al an d the variance to be 
expected due to the imperfect cancellation of the two polarities. Using equations JB It . iB2\ . JB4L and iB5\ . 

E{V} = V Bfj,(B) [E{N + (B)}-E{N-(B)}] = 0, (B6) 

whereas, 

V{V} = [V Bv(B)] 2 (V{N + (B)} + V{N-(B)})=[VoB f i(B)L/l] 2 P(B). (B7) 
The typical Stokes V signals to be expected would be of the order of standard deviation of the distribution of signals ay, 



a v = ^V{V} = (VoL/l) B (j.(B) y/Wfi). (B8) 



5 The linear polarization signals scale with the square of the magnetic field strength and, therefore, the bias toward large field strengths discussed in this 
appendix is even more severe in linear polarization. 
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FIG. B9. — Typical Zeeman signals to be expected for each field strength B. Despite the fact that the weak fields overwhelmingly dominate the distribution 
of field strengths, they leave very small signals. This bias is partly due to the scaling of the Stokes V with the field strength, but it is also due to a preferential 
canceling of the mixed polarity weak fields. The ordinates have been scaled so that the maximum equals one. 



Note that the Zeeman polarization signals due to B do not scale with their filling factor P(B) but with ^/P(B). This sca ling favors 
the contribution of the improbable magnetic fields, that is to say, the contribution of the large field strengths. Figure lB9l shows 
ov normalized so that the maximum equals one. It uses the reference P(B) in Figure[7]plus the ansatz, 



m=\tBBh [(B-750G)/300G] + | ~ { \l 2 \ <™ Q % 



(B9) 



although the latt er is secondary since /x(B) only varies from 1/2, for an isotropic distribution of magnetic field, to 1, for vertical 
fields. Equation dB9l parameterizes our theoretical prejudices; the distribution of inclinations should be more isotropic for weak 
fields (which can be easily tilte d by the drag of the gra nula tion flows) and more vertical for the strong fields (for which the 
buoyancy is important; see, e.g.. lSchiissleJl9 86). Figure lB9l shows how the magnetic fields occupying most of the surface (say, 
B < 200 G; Figure[7} leave Stokes V signals smaller than those of the kG fields. 

In short, it is not difficult to explain why the Zeeman signals tend to detect strong fields, despite the fact that they occupy a small 
fraction of the solar surface. It is due to the scaling of the Stokes V signals with the field strength, plus the random cancellation 
of the mixed polarity. The latter preferentially affects the most probable field strengths, i.e., the smallest ones. 
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